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Using only independent high-scale OPE input, we investigate QCD sum rule constraints 
on two currently incompatible versions of the isovector vector spectral function, one ob- 
tained from electroproduction (EM) data, the other from hadronic r decay data. Sum 
rules involving weighted integrals over the spectral function, from threshold to a variable 
upper endpoint so, are employed. It is shown that both the normalization and slope with 
respect to so of the EM spectral integrals disagree with the corresponding OPE expecta- 
tions, while both normalization and slope are in good agreement when hadronic r decay 
data is used instead. These results favor determinations of the leading hadronic vacuum 
polarization contribution to a M obtained using the r decay data, and hence Standard 
Model predictions for a M compatible with the current experimental determination. 
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1. Introduction 

The recent high-precision experimental determination of = (g — 2)„^has led to 
intense interest in the corresponding Standard Model (SM) prediction^! The largest 
of the non-pure-QED SM contributions, that due to the leading order hadronic vac- 
uum polarization, [a^]£^, can be obtained as an appropriately-weighted integral 
over the electromagnetic (EM) spectral function, Pem( s ) ■ The 1=1 component 
of Pem(s) is related to that of its charged current 1=1 partner, measured in 
hadronic r decay, by CVC. Recent attempts to incorporate r decay data into the 
evaluation of \a^\had using the isospin-breaking (IB) corrected version of this re- 
lation ^ have encountered a consistency problem in the dominant (irir final state) 
contribution to [a^ad- Specifically, even after accounting for known sources of 
IB El the 7T7T component of the resulting r-based 1 = 1 vector spectral function 
is incompatible with that obtained from the high-precision CMD2 e + e~ — > ir + ir~ 
data ®, with surprisingly large (~ 10%) discrepancies remaining in the region of 
7T7T invariant masses from ~ 0.85 to ~ 1 GeV™. SM predictions for a M based solely 
on electroproduction (EM) data indicate a ~ 2.5a deviation from the experimen- 
tal result, while those based o n th e IB-corrected hadronic t decay data indicate 
compatibility at the ~ la level 1^1. The discrepancy thus significantly impacts the 
question of whether or not the current experimental determination of a M ^ shows 
evidence for beyond-the-SM contributions. 
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To resolve the above discrepancy, we have investigated sum rule constraints 
on the EM and r-decay-based data sets, working with sum rules of the form 
f^° h dsw(s) p(s) — -^jk f, , q dsw(s) II(s), where Sth is the relevant threshold, 
II(s) is the correlator with spectral function p(s), w(s) is a function analytic in 
the region of the contour, and the OPE representation of II(s) is to be used on the 
RHS. At intermediate scales, weights satisfying w(s = so) = must be employed, 
in order to suppress duality violating contributions from the vicinity of the timelike 
point™ This is most conveniently done by using the variable y = s/sq and weights 
w(y) satisfying w(l) = 0. Integrated OPE contributions of dimension D then scale 
as 1/sq D 2 ^ 2 , allowing reliable self-consistency checks for the absence of neglected 
higher D contributions. 

In the current study, to be specific, we employ the ALEPH data and covariance 
matrix ^ (with updated normalization) for the hadronic r decay version of the 
isovector vector spectral function. For the various exclusive component s of th e EM 
spectral function and their errors, we follow the assessments of Refs. EDS! More 
recent high-precision data is incorporated as discussed in Rcf. 

ED 

In the region of 

the EM-r discrepancy, we focus on the CMD2 irir data. The weights, w(y), and 
scales, so, used in this study have been chosen in such a way that (i) the OPE is 
essentially entirely dominated by its D = contribution, (ii) the convergence of the 
integrated D = series, order-b y-or der in a s , is excellent, and (hi) poorly known 
D = 6 contributions are absent The OPE sides of the various sum rules are 
then determined, up to very small non-perturbative corrections, by the single input 
parameter, a s (Mz), for which we employ an average of high-scale determinations 
independent of the EM and r data being tested LUJ. The values relevant to lower 
scales are obtained by standard four- loop running and matching H^l. 

For illustration, we present below results for three weight cases, w(y) = 1 — y, 
Ws(y), and ws(y), where Wff(y) = 1 — (-^n") y + ^r- The first weight has a zero of 
order 1 at y = 1, the remaining weights zeros of order 2. The weights are all non- 
negative in the integration region, simplifying the interpretation of the su m ru le 
tests. A discussion of other advantages of these weights may be found in Ref. 



2. Results and Discussion 

We find that the r-decay-based spectral integrals are compatible with the high-scale 
OPE input for a wide range of so, but that the EM data yields spectral integrals 
consistently lower than those predicted by the OPE ^5 The s dependence of the 
OPE and corresponding spectral integrals is also in good agreement for the r c ase, 
but in significant disagreement for the EM case^] Figures may be found in Ref. ^5 
Regarding the normalization of the EM spectral integrals: to quantify the nor- 
malization discrepancy, we have used the spectral integral values at the highest 
accessible common scale for the EM and r cases (so = m r) to obtain an effective 
value for a s {Mz)- This is to be compared to the independent high-scale average, 
a s {M z ) = 0.1200 ± 0.0020. The results arc shown in Table 1. While the values im- 
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plied by the EM data are low by only ~ 2a, typical fluctuations in the data which 
would bring them into better agreement with OPE expectations would also serve to 
increase [o/J^- The r data, in contrast, yields a s (Mz) values in good agreement 
with the high-scale average. Somewhat lower values of a s (Mz) are also easily ac- 
commodated since the overall normalization uncertainty for the r data corresponds 
to an uncertainty of ±0.0010 in the fitted a s (Mz)- 



Table 1. Values of a a (Mz) obtained by fitting to the so = m\ experimental EM and t 
spectral integrals with central values for the D = 2, 4 OPE input. 



Weight 


EM or t 


a s {M z ) 


1-2/ 

W3 

we 


EM 
EM 
EM 


U.lJ.d»_ 0035 

n 1152+ 00019 

u.±±oz_ 0021 


1-2/ 

W3 

we 


r 
r 

T 


n -1910+U.UU2V 
u.izio_ 0032 

u.nao_ 0021 

u.izui_ 00 22 



The problem of the sq dependence of the EM spectral integrals is similarly 
quantified in Table 2, for w(y) = 1 — y, we(y). S exp and Sope are the slopes with 
respect to sq of the spectral and OPE integrals, respectively. OPE entries labelled 
"indep" correspond to the high-scale average a s (Mz) input above, those labelled 
"fit" to the fitted values given in Table 1. The uncertainty in Sope is seen to be very 
small. The very weak dependence on a s means that even the change from "indep" 
to "fit" input on the OPE side has only marginal impact, the OPE versus spectral 
integral slope discrepancy being reduced from 2.6 to 2.3 a for w(y) = 1 — y and 
2.5 to 2.2 a for w(y) — wg(y). The source of the slope problem thus lies entirely on 
the data side. The problem can be cured only by a change in the shape of the EM 
spectral distribution, such as would occur if the IB-corrected r data, rather than 
the EM data, represented the correct version of p£m{ s )- 



Table 2. Slopes wrt so of the EM OPE and spectral integrals 



Weight 


Sexp 


a 3 (M z ) 


Sope 


1-2/ 


.00872 ± .00026 


indep 


.00943 ± .00008 






fit 


.00934 ± .00008 


we 


.00762 ± .00017 


indep 


.00811 ± .00009 






fit 


.00805 ± .00009 



We conclude that either (i) non-one-photon physics effects, as yet unidentified, 
are contaminating the EM data or (ii) there remain experimental problems with 
the (pre-2005) EM data. The latter conclusion is favored by the just-released SND 
e + e~ — » 7r + 7r~ results ED which are compatible with the IB-corrected hadronic r 
decay data. In either case it follows that determinations of [a M ]j^ incorporating r 
decay data are favored over those based solely on EM cross-sections and hence that 
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the SM prediction for is in good agreement with the current experimental result. 

A final comment concerns the IB correction associated with p-Lu "mixing" . 
The current standard evaluation El is based on a chirally-constrained model (the 
"GP/CEN model") fitted to a since-corrected form of the CMD2 data. Because 
of strong cancellations, the model dependence of the integrated [a^ad interference 
contributio n tu rns out to be significantly larger than any one model's fitting-induced 
uncertainty^! Table 3 shows the results of fits to the CMD2 data for the Gounaris- 
Sakurai (GS), hidden local symmetry (HLS), and Kuhn-Santamaria (KS) models, 
as well as two versions of the GP / CEN model one a refitting to the current 
version of the CMD2 data (GP/CEN^), the other a refitting which incorporates an 
additional phase, such as arises when direct IB lu — > tttt decay contributions are 
taken into account (GP/CEN*). Smaller values of ^(a^)] w also serve to reduce 
the difference between the EM and r-based determinations of [dju] fi- 



nable 3. Model dependence of the integrated p-uj interference IB correction 



Model 


X 2 /dof 


[SM] p _ u x 10"' 


GS 


35.9/38 


2.0 ±0.5 


HLS 


36.6/38 


4.0 ±0.6 


KS 


37.1/38 


3.8 ±0.6 


GP/CEN* 


40.6/39 


2.0 ±0.5 


GP/CENt 


61.5/40 


3.7 ±0.7 
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